arXiv: 1501.00163vl [cs.IT] 31 Dec 2014 


l 


Error Correction in Polynomial Remainder Codes 
with Non-Pairwise Coprime Moduli and Robust 
Chinese Remainder Theorem for Polynomials 

Li Xiao and Xiang-Gen Xia 


Abstract —This paper investigates polynomial remainder codes 
with non-pairwise coprime moduli. We first consider a robust 
reconstruction problem for polynomials from erroneous residues 
when the degrees of all residue errors are assumed small, namely 
robust Chinese Remainder Theorem (CRT) for polynomials. It 
basically says that a polynomial can be reconstructed from 
erroneous residues such that the degree of the reconstruction 
error is upper bounded by r whenever the degrees of all residue 
errors are upper bounded by r, where a sufficient condition for 
t and a reconstruction algorithm are obtained. By releasing the 
constraint that all residue errors have small degrees, another 
robust reconstruction is then presented when there are multiple 
unrestricted errors and an arbitrary number of errors with small 
degrees in the residues. By making full use of redundancy in 
moduli, we obtain a stronger residue error correction capability 
in the sense that apart from the number of errors that can be 
corrected in the previous existing result, some errors with small 
degrees can be also corrected in the residues. With this newly 
obtained result, improvements in uncorrected error probability 
and burst error correction capability in a data transmission are 
illustrated. 

Index Terms —Burst error correction, error correction codes, 
polynomial remainder codes, residue codes, robust Chinese Re¬ 
mainder Theorem (CRT). 

I. Introduction 

HE Chinese Remainder Theorem (CRT) can uniquely 
determine a large integer from its remainders with respect 
to several moduli if the large integer is less than the least 
common multiple (1cm) of all the moduli JT|, (2). Based on the 
CRT for integers, residue codes with pairwise or non-pairwise 
coprime moduli are independently constructed, where code¬ 
words are residue vectors of integers in a certain range modulo 
the moduli. More specifically, a residue code with pairwise co¬ 
prime moduli mi, • ■ ■ , rrik, rrik+i, • ■ ■ , m n consists of residue 
vectors of integers in the range [0, n^=i TO i)> where the first k 
moduli form a set of nonredundant moduli, and the last n — k 
moduli form a set of redundant moduli used for residue error 
detection and correction. Over the past few decades, there has 
been a vast amount of research on residue error correction 
algorithms for such a class of codes. For more details, we refer 
the reader to 0-0. By removing the requirement that the 
moduli be pairwise coprime, a residue code with non-pairwise 
coprime moduli ■ ■ ■ ,mi consists of residue vectors 

of integers in the range [0, lcm(mi, m 2 , • • • , mi)). Compared 
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with residue codes with pairwise coprime moduli, the residue 
error detection and correction algorithm for residue codes with 
non-pairwise coprime moduli is much simpler, and the price 
paid for that is an increase in redundancy. Moreover, residue 
codes with non-pairwise coprime moduli may be quite effec¬ 
tive in providing a wild coverage of “random” errors p24)— [26]. 
In order to perform reliably polynomial-type operations (e.g., 
cyclic convolution, correlation, DFT and FFT computations) 
with reduced complexity in digital signal processing systems, 
residue codes over polynomials (called polynomial remainder 
codes in this paper) with pairwise or non-pairwise coprime 
polynomial moduli have been investigated as well [291-1361. 
where codewords are residue vectors of polynomials with 
degrees in a certain range modulo the moduli and all poly¬ 
nomials are defined over a Galois field. Polynomial remainder 
codes are a large class of codes that include BCH codes 
and Reed-Solomon codes as special cases |27|, (281. Due to 
two important features in residue codes: carry-free arithmetics 
and absence of ordered significance among the residues, 
residue error detection and correction technique in residue 
codes has various applications in, for example, fault-tolerant 
execution of arithmetic operations in digital processors and in 
general digital hardware implementations on computers 01- 
El, ED, ED, orthogonal frequency division multiplexing 
(OFDM) and code division multiple access (CDMA) based 
communication systems un-urn and secure distributed data 
storage for wireless networks 1201 - 1231 . 

In this paper, we focus on polynomial remainder codes with 
non-pairwise coprime moduli. Note that a coding theoretic 
framework for such a class of codes has been proposed in (33}, 
where the concepts of Hamming weight, Hamming distance, 
code distance in polynomial remainder codes are introduced. 
It is stated in [33|] that a polynomial remainder code with non- 
pairwise coprime moduli and code distance d, can correct up 
to [(d — 1)/2J errors in the residues, and a fast residue error 
correction algorithm is also presented, where |* is the floor 
function. This reconstruction from the error correction method 
is accurate but only a few of residues are allowed to have errors 
and most of residues have to be error-free. The goal of this 
paper is to study robust reconstruction and error correction 
when a few residues have arbitrary errors (called unrestricted 
errors) similar to l33l and some (or all) of the remaining 
residues have small errors (i.e., the degrees of errors are small). 
It is two-fold. One is to study robust reconstruction and the 
other is to study error correction, i.e., accurate reconstruction, 
when residues have errors. 
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Considering instabilities of data processing in wireless sen¬ 
sor networks and signal processing systems, robust reconstruc¬ 
tions based on the CRT for integers were recently studied 
in El-® with different approaches. In this paper, by 
following the method in El together with the error correction 
algorithm for polynomial remainder codes in El . we first 
propose a robust reconstruction algorithm for polynomials 
from erroneous residues, called robust CRT for polynomials, 
i.e., a polynomial can be reconstructed from erroneous residues 
such that the degree of the reconstruction error is upper 
bounded by the robustness bound r whenever the degrees of 
all residue errors are upper bounded by r, where a sufficient 
condition for r for the robustness to hold is obtained. Next, 
by releasing the constraint that the degrees of all residue 
errors have to be bounded by r, we propose another robust 
reconstruction algorithm when a combined occurrence of 
multiple unrestricted errors and an arbitrary number of errors 
with degrees upper bounded by A happens to the residues, 
where a sufficient condition for A is also presented in this 
paper. Note that a combined occurrence of a single unrestricted 
error and an arbitrary number of small errors in the residues 
was considered for the robust reconstruction based on the CRT 
for integers in El . but its approach is hard to deal with the 
case of multiple unrestricted errors combined with small errors 
in the residues due to a considerable decoding complexity. A 
detailed comparison in terms of robust reconstruction between 
this paper and El . f39l is pointed out later in this paper (see 
Remark [3J. One can see that the above reconstructions may 
not be accurate but robust to the residue errors in terms of 
degree and all the residues are allowed to have errors. 

Finally, we consider the residue error correction in a poly¬ 
nomial remainder code with code distance d. Compared with 
the result in El . by making full use of the redundancy in 
moduli and newly proposed robust reconstruction method, we 
obtain a stronger residue error correction capability in the 
sense that apart from correcting up to |_(d — 1)/2J residue 
errors, a polynomial remainder code with code distance d can 
correct some additional residue errors with small degrees. With 
this newly obtained result, improvements in the performances 
of uncorrected error probability and burst error correction 
considered in a data transmission are illustrated. 

The rest of the paper is organized as follows. In Section 
m we briefly introduce some fundamental knowledge in poly¬ 
nomials over a Galois field and coding theory of polynomial 
remainder codes with non-pairwise coprime moduli obtained 
in El In SectionlnU we propose robust CRT for polynomials. 
In Section [IV] another robust reconstruction is considered 
when a combined occurrence of multiple unrestricted errors 
and an arbitrary number of errors with small degrees is in 
the residues. In Section [V] a stronger residue error correction 
capability in polynomial remainder codes with non-pairwise 
coprime moduli and its improvements in uncorrected error 
probability and burst error correction in a data transmission 
are presented. We conclude this paper in Section [VI] 

II. Preliminaries 

Let F be a field and F[x\ denote the set of all polynomials 
with coefficients in F and indeterminate x. The highest power 


of x in a polynomial fix) is termed the degree of the 
polynomial, and denoted by deg (f(x)). All the elements of F 
can be expressed as polynomials of degree 0 and are termed 
scalars. A polynomial of degree n is called monic if the co¬ 
efficient of x n is 1. Denote by gcd (j\ (x), f' 2 (x). ■ ■ ■ , /l(x)) 
the greatest common divisor (gcd) of a set of polynomials 
fi(x), i.e., the polynomial with the largest degree that divides 
all of the polynomials fi(x) for 1 < i < L. The least 
common multiple (1cm) of a set of polynomials /,(x), denoted 
by 1cm {fi{x),f 2 (x), • • • , /l(x)), is the polynomial with the 
smallest degree that is divisible by every polynomial fi(x) for 
1 < i < L. For the uniqueness, gcd(-) and lcm(-) are both 
taken to be monic polynomials. Two polynomials are said to 
be coprime if their gcd is 1 or any nonzero scalar in F. A 
polynomial is said to be irreducible if it has only a scalar 
and itself as its factors. The residue of f(x) modulo g(x) is 
denoted as [f(x)] g r x y Throughout the paper, all polynomials 
are defined over a field F , and *j and [*] are well known 
as the floor and ceiling functions. 

Let mi(x),m2(x), ■ ■ ■ , rriL{x) be L non-pairwise coprime 
polynomial moduli, and M(x ) be the 1 cm of all the mod¬ 
uli, i.e., M(x) = lcm(mi(a;),7712(2),• ■ • , TOl(x)). For any 
polynomial a(x) with deg(a(x)) < deg(M(x)), it can be 
represented by its residue vector (ai(x), as2(x), • • ■ ,cll{x)), 
where a,(a) = [a(x)] m . (a;) , i.e., 

a(x) = ki(x)rrii( x) + a*(x) (1) 


with deg (a,(x)) < deg (nii(x)) and ki(x) £ F[x] for 1 < * < 

L. Here, we call such k, (x) in ([TJ the folding polynomials. 
Equivalently, a(x) can be computed from its residue vector 
via the CRT for polynomials 0], El . 


a(x) 


" L 

Y, Oi(x)Di(x)Mi(x) 

_i=1 


M(x) 


( 2 ) 


where M,;(x) = , _D,(x) is the multiplicative inverse of 

Mi(x) modulo Hi(x), if /^(x) ^ 1, else D^x) = 0, and 
{/ 2 i(x)}((_ 1 is a set of L pairwise coprime monic polynomials 
such that nf =1 /L(x) = M(x) and /Xj(x) divides nii(x) for 
each 1 < i < L. Note that if n~ii(x) are pairwise coprime, 
we have //j(x) = m;(x) for 1 < i < L, and then the above 
reconstruction reduces to the traditional CRT for polynomials. 

As seen in the above, polynomials a(x) with deg (a(x)) < 
deg (M(x)) and their residue vectors are isomorphic. Further¬ 
more, the isomorphism holds for the addition, substraction, 
and multiplication between two polynomials a(x) and b(x), 
both with degrees less than deg (M(x)). First convert each 
polynomial to a residue vector as 


a(x) <-)• (ai(x), • • • ,a L {x)) and b{x) -H- (6i(x),--- ,b L (x)). 

(3) 

Then, the residue representation of c(x) = a(x) ± b(x) or 
d(x) = \a{x)b{x)]M( x ) is given by, respectively, 

c(x) (ai(x) ± bi(x), • • • , ol(x) ± b L (x)) , (4) 

d{x) fA ([ai(x)6i(x)] mi(x) , • • • , [a L (x)b L (x)] mL ( x) ) . (5) 

Moreover, an important property in a polynomial remainder 
code with non-pairwise coprime moduli is that if a(x) = 
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a,i(x) mod nii(x) and a(x) = dj(x) mod rrij(x), the follow¬ 
ing congruence holds (T), ||33l : 

a,i{x) = a,j(x) mod dij(x), (6) 

where dij(x) = gcd (wj(ir), m,j(x)). We call equation <[6]» 
a consistency check between residues at(x) and a,j(x). If 
<[6]» holds, a,i(x) is said to be consistent with aj(x); other¬ 
wise, di{x) and cij(x) appear in a failed consistency check. 
A residue vector (ai(ir), 02 ( 2 ), • ■ ■ ,cil(x)) is said to be a 
polynomial remainder codeword if it satisfies the consistency 
checks given by © for all pairs of residues in the vector. 
So, any polynomial a(x ) with deg(a(cc)) < deg(M(x)) is 
represented by a unique polynomial remainder codeword, i.e., 
its residue vector. Conversely, every polynomial remainder 
codeword is the representation of a unique polynomial with de¬ 
gree less than deg (M(x)). We call the set of such codewords 
a polynomial remainder code with moduli m \ (x), • • • , rrir (x), 
which is linear according to (0]). 

If t errors, (x), • • • ,ei t (x), in the residues have occurred 
in the transmission, then the received residue vector, denoted 
by (di(ir), • • • ,o,l{x)), is determined by 

(ai(x), ■ • ■ , a L (x )) = (ai(x), • ■ • , a L (x)) (7) 

+ (0, • • • , e h (x), • • • , e i2 (x), ■ ■ ■ ,e it (x),---), 

where deg (e,. (x)j < deg (mu (s)) for 1 < j < t, and 
the subscripts «],••• ,i t are the corresponding positions of 
the residue errors (*),-•• ,ei t (x). In (33) . the capability 
of residue error correction in a polynomial remainder code 
with non-pairwise coprime moduli has been investigated, and a 
simple method for residue error correction has been proposed 
as well. Before briefly reviewing them, let us present some 
notations and terminologies in polynomial remainder codes 
with non-pairwise coprime moduli used in (33). Hamming 
weight of a codeword is the number of nonzero residues in 
the codeword, Hamming distance between two codewords is 
defined as the Hamming weight of the difference of the two 
codewords, and code distance of a polynomial remainder code 
is the minimum of the Hamming distances between all pairs of 
different codewords. Due to its linearity ©), the code distance 
is actually equal to the smallest Hamming weight over all 
nonzero codewords. Similar to a conventional binary linear 
code, a polynomial remainder code with code distance d can 
detect up to d - 1 errors in the residues, and correct up to 
|_(d — 1)/2J errors of arbitrary values in the residues. A test 
for the code distance of a polynomial remainder code with 
non-pairwise coprime moduli is presented in the following. 

Proposition 1: (33) Let rrii(x), 1 < i < L, be L non- 

pairwise coprime polynomial moduli, and denote M(x) = 
lcm (mi (x), rri 2 (x), •• • ,rriL(x)). Write M{x) in the form 

M(x) = pi(x) tl p 2 {x) t2 ■ ■■p K {x) tK , (8) 

where the polynomials Pi(x) are pairwise coprime, monic and 
irreducible, and t, is a positive integer for all 1 < i < K. 
For each 1 < i < K, let di represent the number of 
moduli that contain the factor p^x)^. Then, the code distance 
of the polynomial remainder code with the set of moduli 
m 2 (x), ■■ ■ , rriL^x )} is d = min{di, d 2 , • • • , dn}- 


Based on Proposition Q] an explicit method of constructing 
a polynomial remainder code with code distance d is also 
proposed in (33). Let M(x) be decomposed into the product 
of several smaller, pairwise coprime, and monic polynomials 
Pi( x) ti as in the form ©, L represent the number of moduli in 
the code, and d be a positive integer such that 1 < d < L. For 
each pi{ a:)* 4 , assign pi(x) ti to di different moduli, such that 
di > d, with the equality for at least one i. Set each modulus 
to be the product of all polynomials assigned to it. Then, 
the resulting polynomial remainder code will have the code 
distance d. In particular, repetition codes can be obtained in 
the above construction by setting di = d = L, i.e., all moduli 
are identical. Next, the polynomial remainder code defined in 
Proposition © enables fast error correction, as described in the 
following propositions. 

Proposition 2: (33) In a polynomial remainder code with 
code distance d defined in Proposition Q] if only t < L(d- 

1)/2J errors in the residues have occurred in the transmission, 
each erroneous residue will appear in at least [(d — l)/2] +1 
failed consistency checks. In addition, each correct residue will 
appear in at most [(d — 1)/2J failed consistency checks. 

Proposition 3: (33) Let moduli m.i{x) for 1 < 1 < L, 

M(x) and d be defined in Proposition Q] Then, the least 
common multiple of any L — (d—l) moduli is equal to M(x). 

Based on Propositions [2] 0 a polynomial remainder code 
with code distance d can correct up to [(d — 1)/2J residues 
errors, i.e., a(x) can be accurately reconstructed from all the 
error-free residues that can be fast located through consistency 
checks for all pairs of residues di(x) for 1 < i < L. With 
the above result, it is not hard to see the following decoding 
algorithm for polynomial remainder codes with non-pairwise 
coprime moduli and code distance d. 

1) Perform the consistency checks by © for all pairs of 
residues ai(x), 1 < i < L, in the received residue vector. 

2) Take all of those residues each of which appears in at 
most [(d— 1)/2J failed consistency checks. If the number 
of such residues is zero, i.e., for every i with 1 < i < L, 
&i(x) appears in at least [(d—1)/2J +1 failed consistency 
checks, the decoding algorithm fails. Otherwise, go to 3). 

3) If all the residues found in 2) are consistent with each 
other, use them to reconstruct a(x) as a(x) via the 
CRT for polynomials in (2). Otherwise, a(x) cannot be 
reconstructed and the decoding algorithm fails. 

According to Propositions©© if there are [(d — 1)/2J or 
fewer errors in the residues, a(x) can be accurately recon¬ 
structed with the above decoding algorithm, i.e., a(x) = a(x). 
However, if more than [(d— 1)/2J errors have occurred in the 
residues, the decoding algorithm may fail, i.e., a(x) may not 
be reconstructed, or even though a(x) can be reconstructed 
as a(x), a(x ) = a(x) may not hold. In the rest of the paper, 
we assume without loss of generality that the non-pairwise 
coprime moduli toi(ie), m 2 (x), ■ ■ ■ ,itil{x) are L arbitrarily 
monic and distinct polynomials with degrees greater than 0, 
and the following notations are introduced for simplicity: 

1) dij(x) = gcd (mi(x),m,j(x)) for 1 < i,j < L,i ± j\ 

2) Tij = deg (dij(x)) for 1 < i,j < 

3) Tj = min{Tjj, for 1 < * < L,i ^ j}\ 

l 
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4 ) r «W = for 1 <i,j<L,i^ j\ 

5) v/' l! denotes the code distance of the polynomial remain¬ 
der code with moduli T,,i(x) for 1 < j < L, j ^ i, which 
can be calculated according to Proposition [Q 

6 ) ri(i) denotes the v'-th smallest element in an array of 
positive integers S = {ni, n 2 , ■ ■ ■ , uk}- It is obvious that 
n (i) = miniS and n^K) = max S. See S = {3,1, 2, 3, 8 } 
for example, and we have = 1 , ri( 2 ) = 2 , 77 .( 3 ) = 

71 ( 4 ) = 3 , 71 ( 5 ) = 8 . 


III. Robust CRT for polynomials 

Let rriiix), 1 < i < L, be L non-pairwise coprime 
polynomial moduli, M (x) be the 1cm of the moduli, and d be 
the code distance of the polynomial remainder code with the 
moduli. As stated in the previous section, a polynomial a(x ) 
with deg ( a(x )) < deg ( M(x )) can be accurately reconstructed 
from its erroneous residue vector (ai(x), • • • , cll(x)), if there 
are [( d — 1)/2J or fewer errors affecting the residue vector 
(ai(x),--- , az,(x)). Note that the reconstruction of a(x) is 
accurate but only a few of the residues are allowed to have 
errors, and most of the residues have to be error-free. In this 
section, we consider a robust reconstruction problem on which 
all residues a* * (a;) for 1 < i < L are allowed to have errors 
ei(x) with small degrees. 

Definition 1 ( Robust CRT for Polynomials): Q A CRT for 
polynomials is said to be robust with the robustness bound r 
if a reconstruction a(x) can be calculated from the erroneous 
residues cq(x) for 1 < i < L such that deg(a(x) — a(x)) < r 
whenever the residues are affected by errors with degrees 
upper bounded by r, i.e., deg(ej(x)) < r < deg (mfix)) 
for 1 < i < L. 

This robust reconstruction problem we are interested in is 
two-fold: one is how we can robustly reconstruct a(x); the 
other is how large the robustness bound r can be for the 
robustness to hold. The basic idea for the robust CRT for 
polynomials is to accurately determine one of the folding 
polynomials. Consider an arbitrary index j with 1 <j<L. 
If the folding polynomial kj(x) is accurately determined, a 
robust estimate of a(x) can then be given by 

a(x) = kj (x)rrij (x) + dj (x) 

= kj(x)rrij(x) + aj(x) + ej(x), 

i.e., deg (a(x) — a(x)) = deg (e.j(x)) < t. Therefore, the 
problem is to derive conditions under which kj(x ) can be 
accurately determined from the erroneous residues a, (x) for 
1 < i < L. To do so, we follow the algorithm in [!39l| for 
integers. 

Without loss of generality, we arbitrarily select the first 
equation or remainder for i = 1 in © as a reference to 
be subtracted from the other equations for 2 < i < L, 

1 The general robustness is that the reconstruction error is linearly bounded 
by the error bound r of the observation. It is well known that the traditional 
CRT (with pairwise coprime moduli) is not robust in the sense that a small 
error in a remainder may cause a large reconstruction error m, ej. 


respectively, and we have 

{ ki{x)m\{x) — fc 2 (x)TO 2 {x) = a 2 (x) — ai(x) 
k 1 (x)mi(x) - k 3 (x)m 3 (x) = a 3 (x) - a 3 (x) 

. 00 ) 

k 1 (x)m 1 (x) - k L {x)m L (x) = a L {x) - ai(x). 

Denote qn(x) = for 2 < i < L. Then, dividing 

du(x) from both sides of the (i — l)-th equation in dTOb for 
2 < 7 < L, we can equivalently write © as 

ki(x)r 12 (x) - k 2 (x)r 2 i(x) = q 2 i(x) 
ki(x)r 13 (x) - fc 3 (x)r 3 i(x) = g 3 i(x) 

: (11) 
h(x)T 1L (x) - k L (x)V L1 (x) = q L i(x). 

Since 1 1 7 ; (rr;) and I ’, | (x) are coprime, by Bezout’s lemma for 
polynomials we have 


k\{x) = g J ;i(x)T li (x) + fc(x)g i i(x)r il (x), for 2 < i < L, 

( 12 ) 

where k(x) is some polynomial in F[x], and Tij(x) is 
the multiplicative inverse of Tii(x) modulo Tii(x), i.e., 
f H (x)r H (x) = 1 mod Tji (x). 

Next, we can use 


q% i (x) = 


afix) - di(x) - [di(x) - ai(x)] di . 


(x) 


= qn(x)+ 


du(x) 

e i(x) - ei(x) - [efix] - ei(x)] dii(a .) 


du(x) 


03) 


as an estimate of qn(x) for 2 < i < L in (fl2ll . and we have 
the following algorithm. 


Algorithm I: 


Step 1: Calculate du(x) = gcd (mi(x), mfix)). 


r H (x) = 


mi ( x ) 
dn(x) 


and Tji(x) = 


(a) 

du(x ) 


for 2 < i < L 


from the given moduli rrij (x) for 1 <3< L, which can 
be done in advance. 

Step 2: Calculate qn(x) for 2 < i < L in (IT3l > from 
m,j(x) and the erroneous residues aj(x) for 1 < j < L. 
Step 3: Calculate the remainders of < 7 ii(x)Tii(x) modulo 
Tii(x), i.e.. 


^ii(x) = qn{x)T u (x) mod Tji(x) 


(14) 


for 2 < i < L, where f’],;(x) is the multiplicative 
inverse of Tii(x) modulo T,i(x) and can be calculated 
in advance. 

• Step 4: Calculate fci(x) from the following system of 
congruences: 

fei(x) = (x) mod Ta(x), for 2 < i < L, (15) 


where moduli Ta (x) may not be pairwise coprime. Note 
that fci(x) is calculated by using the decoding algorithm 
for the polynomial remainder code with moduli Ta(x) 
for 2 < i < L, based on Propositions [3 0 in Section nn 
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Remark 1: As we mentioned before, the basic idea in the 
robust CRT for polynomials is to accurately determine one 
of folding polynomials, which is different from the robust 
CRT for integers l38l . |[39l where all folding integers are 
accurately determined and each determined folding integer 
provides a reconstruction, and all the reconstructions from 
all the determined folding integers can then be averaged to 
provide a better estimate. Accordingly, since we do not need 
to calculate other folding polynomials ki(x) for 2 < i < L in 
the above Algorithm /, qn{x) = qn(x) in (fl3l > does not have 
to hold for all 2 < i < L, that is, residues £n(x) in (IT5l) . 
2 < i < L, are allowed to have a few errors. This is why 
we use the decoding algorithm in Section [IT] to reconstruct 
k\(x) in Step 4 in terms of the polynomial remainder code 
with moduli IVi (x) for 2 < i < L. 

Let w/ 1 ) denote the code distance of the polynomial re¬ 
mainder code with moduli I Yi (x) for 2 < i < L and r be the 
robustness bound, i.e., deg (ej(x)) < r for 1 < i < L. With 
the above algorithm, we have the following lemma. 

Lemma 1: k\{x) can be accurately determined in Algorithm 
I, i.e., k\(x ) = ki(x), if the robustness bound r satisfies 

T < T 1(L(^ (1) -1)/2J+1)> (16) 

where = deg (d 1? (a;)) for 1 < i, j < L, i ^ j , and r 1(j ) 
denotes the j-th smallest element in {ti 2 ,ti 3 , • • • , til}. 

Proof: Let T(a;) = lcm (T 2 i(x), T 3 i(x), • • • , T L i(x)). 
We first prove deg ( ki(x )) < deg (T(x)). If deg ( ki(x )) = 0, it 
is obvious for deg (, ki(x )) < deg (r(x)). If deg (ki(x)) ^ 0, 
we have deg {k\(x)m\(x)) = deg(fci(x)) + deg (mi (a;)) = 
deg (a(x)) < deg (M(x)) = deg (mi(a;)r(x)) = 

deg (mi(a;))+deg (T(a;)), and thus deg (ki(x)) < deg (T(a;)). 

Among the residue errors Bi(x) for 2 < i < L, there are 
v > L — 1 — [ u)l ,^~ 1 j errors for 1 < j < v such 

that deg (e^. (x) — ei(a;)) < r < deg (du^x)) according to 
©• So, we have [ei } (x) - ei(x)\ dli . (x ) = e^(x) - ei(x). 
From ©, we have qi j 1 ( 1 ) = qi.i(x), and it is not hard 
to see from (ThB that k\(x) and fa.\ (x)Tij.(a;) have the same 
remainder modulo Tj -i(a;), i.e., k\{x) = £j.i(x) mod Li i(a;). 
Then, regard (£21 (a:), £31 ( 2 :), • • • ,£ 7 , 1 ( 2 ;)) in ([15]) as an er¬ 
roneous residue vector of k\{x) modulo Tji(s) for 2 < 
i < L. Since there are at most [(i// 1 ) — 1)/2J errors in 
(£2i(y)j £ 31 ( 2 ;), • • • , £li(*)), we can accurately determine the 
folding polynomial k\{x) in Step 4 of Algorithm /, i.e., 
k\(x) = k\(x), by applying the residue error correction 
algorithm based on Propositions [2j El for the polynomial 
remainder code with moduli T;i(:r) for 2 < i < L. ■ 

Recall that m\(x) or ai(a;) in the above Algorithm I is 
arbitrarily selected to be a reference, which is not necessary. 
In fact, any remainder can be taken as the reference. In 
order to improve the maximal possible robustness bound, we 
next present the following theorem through selecting a proper 
reference folding polynomial. 

Theorem 1: If the robustness bound r satisfies 

T < i<£c X l { r i(L(™«-l)/2J+l)} . ( 1? ) 

where w^ 1 ’ is the code distance of the polynomial remainder 
code with moduli T'ji(x) for 1 < j < L,j i, and ruj) 


denotes the j-th smallest element in {r^, for 1 < k < L, k 
i}, then a(x) can be robustly reconstructed through Algorithm 
L that is, the robust CRT for polynomials in Definition[T]holds. 

Proof: Let us choose such an index i(, that 

T i 0 (L(i« (i t |) -l)/ 2 j+l) = { T i(L(^ (i) -l)/ 2 j+l)} ’ (I**) 

Then, replacing the index 1 with i$ and taking d-, 0 (x) as the 
reference in Algorithm /, we can accurately determine k, 0 (x) 
under the condition ( fTTI ). thereby robustly reconstructing a(x) 
as a(x) in (©, i.e., deg (a(x) — a{x)) < r. ■ 

Remark 2: From ( fl 8 ] i. it guarantees that there are at most 
L(ti)(®°) — 1 )/ 2 J errors in the residues £,u 0 {x) for 1 < i < 
L,i io in Step 4 of Algorithm I to determine ki 0 {x) with 
respect to moduli ra 0 (x). If > 3, ki 0 {x) is accurately 
determined based on the residue error correction algorithm 
in Section m for the polynomial remainder code with moduli 
Fjj 0 ( 2 ;) for 1 < i < L, i i$. If w^ l °' > < 3, all the residues 
fa 0 ( x ) f° r k < i < L,i io, are accurate and consistent, and 
ki 0 (x) is accurately determined via the CRT for polynomials 
from all the L — 1 residues £ii 0 (x). 

Example 1: Let us consider a notable class of polyno¬ 
mial remainder codes with special moduli (the corresponding 
integer residue codes were introduced in l24ll . l37l ). i.e., 
m,i(x ) = ]”[ dij(x) for 1 < i < L and {dij(x), for 1 < 

fe[i 

i < L\i < j < L} are pairwise coprime. Let ^ 12 ( 2 ;) = 
(x + l) 4 , d 13 (x) = (x - l) 4 , du(x) = (x + 2) 4 , d 23 (x) = 
(x — 2) 4 ,d24{x) = (x + 3) 4 , d 3 i(x) = {x — 3) 4 . Since 
deg ( dij(x )) = 4 holds for every 1 < i, j < 4, i j, it is easy 
to see from Theorem [I] that the robustness bound is r < 4, i.e., 
any a(a;) with deg (a(x)) < deg (lcm(mi(x), ■ ■ ■ , 7714 ( 2 ;))) = 
24 can be robustly reconstructed from its erroneous residues 
when the degrees of all residue errors are less than 4. 

If the above result is referred to as the single stage robust 
CRT for polynomials, multi-stage robust CRT for polynomials 
can be easily derived by following the method used for integers 
in {39) • Similarly, multi-stage robust CRT for polynomials may 
improve the bound for r obtained in Theorem [T] for a given 
set of polynomial moduli. Another remark we make here is 
that a residue error ei(x ) is said to be a bounded error with 
an error bound l if its degree is less than or equal to /, where 
l is a small positive integer. What Theorem Q] tells us is that 
for the set of moduli Wj( x) in the above, a polynomial a(x) 
with deg(a(o;)) < deg(M(x)) can be robustly reconstructed 
from its erroneous residues if all residue errors are bounded, 
and the error bound r is given by dm Later, the constraint 
that all residue errors are bounded will be released, and the 
combined occurrence of multiple unrestricted errors and an 
arbitrary number of bounded errors in the residues will be 
considered in the next section. 

IV. Robust reconstruction under multiple 

UNRESTRICTED ERRORS AND AN ARBITRARY NUMBER OF 
BOUNDED ERRORS IN THE RESIDUES 

Consider again the L non-pairwise coprime moduli rrii(x) 
for 1 < i < L. In this section, we assume that there are 
t < L^ - 1 )/ 2 J unrestricted errors and an arbitrary number of 
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bounded errors with the error bound A in the received residue 
vector (di(x),--- , dz,(x)). Similarly in this case, the robust 
reconstruction problems for us are: 1 ) how can we robustly 
reconstruct a(x)7 2) how large can the error bound A be for 
the robustness to hold? Note that d > 3 is necessarily assumed 
in this section, otherwise it is degenerated to the case of robust 
CRT for polynomials in Section [III] Therefore, due to the 
existence of unrestricted residue errors, the bound for A is 
expected to be smaller than or equal to the bound for r as in 
(ITTI) . In order to answer the above questions, we first give the 
following lemmas. 

Lemma 2: Let ur- ' > denote the code distance of the polyno¬ 
mial remainder code with moduli T(x) for 1 < j < L, j ^ i. 
Then, we have 

minjw/ 1 ), w^ 2 \ ■■■ , w = d, (19) 

where d is the code distance of the polynomial remainder code 
with moduli TOj(x) for 1 < i < L. 

Proof: First, let us prove w' 1 ' 1 > d for each 1 < i < L. 
Without loss of generality, we only need to prove vr 1 ' 1 > d. 
Let M(x) be written as in (|S), i.e., 

M(x) = pi(x) tl p 2 {x) t2 • • ■p K {x) tK , (20) 

where the polynomials Pi(x) are pairwise coprime, monic 
and irreducible, and t, is a positive integer for all 1 < % < 
K. Define F(x) = lcm(r 2 i(x),F 3 i(x), • • • ,Fii(x)). Since 
M(x) = mi(x)r(x), we can write TOi(x) and T(x) as 

to \{x) = pi(x) ll p 2 (x) 12 ■ ■ ■ pk{x) 1k and 

T(x) = pi(x) tl_/l p 2 (x )* 2_Z2 ■ ■ -p K {x) tK ^ lK , ( 21 ) 

where U > k > 0 for 1 < i < K. First, consider 
0 < U < U, and let w \ 1 ' represent the number of moduli 
r 2 i(x), • • • ,Fii(a;) that contain the factor pi{ x) ti ~ li . In this 
case, we have u:[ 11 = di, where d, is defined in Proposition Q] 
because for every j with 2 < j < L, rrij{x) contains Pi{x) ti 
if and only if its corresponding Fji(x) = contains the 

factor Next, consider /,■ = f,, and r(x) does not 

contain the item of Pi[x). Hence, according to Proposition [T] 
is the minimum of {di,d 2 , • • • ,d k}. i.e.. = d, if 

0 < U < ti holds for all 1 < i < K, else is the minimum 
of a subset of { di , d 2 , • ■ ■ , cZr-}, i.e., > d. So, we have 

W' 11 > d. Note that the above proof is independent of an 
arbitrary choice i = 1 for u;W. Therefore, we have ti>M > d 
for each 1 < i < L. 

Next, we prove that there is at least one i such that = d. 
Without loss of generality, we assume that d\ = d. From the 
above analysis, if if/ 1 ) > d , we must have i\ = t\, i.e., mi(x) 
contains the factor pi(x) tl . Similarly, if all for 1 < i < L 
are strictly larger than d, we know that all m\ (x) for 1 < i < L 
contain the factor p^xf 1 , i.e., d\ = d = L. Thus, di = L 
for all 1 < i < L. It is in contradiction with the assumption 
in the end of Section El that toi (x), m 2 (x), ■ • ■ ,m,L(x ) are 
monic and distinct polynomials with degrees greater than 0. 
Thus, we have minju/ 1 ), • • • , w^} = d. ■ 

Lemma 3: Let d denote the code distance of the polynomial 
remainder code with moduli m.i(x) for 1 < i < L. Assume 
that there are t < ~ 1 )/ 2 J unrestricted errors, and 


any other error is bounded in the received residue vector 
(ai(x), ■ ■ ■ ,o,l(x)). The error bound A here is assumed less 

than Ti, where t\ = min{r 1:) ', for 2 < j < L} and 

j 

t ij = deg (dij(x)). If di(x) is known as an error-free residue 
or a residue with a bounded error, we can accurately determine 
k\{x) using Algorithm L i.e., fci(x) = k\(x). However, if 
di(x) is known as a residue with an error of degree greater 
than A, k\ (x) may not be reconstructed, and even though k\ (x) 
is reconstructed in Algorithm /, k\{x) = fci(x) may not hold. 

Proof: If u n (x) with i > 1 is an error-free residue or 
a residue with a bounded error, i.e., e^(x) = 0 or ej(x) 

0 with deg(e,(x)) < A, we have e^x) — e±(x) = [e.j (x) — 
ei(x)] rf i_( x ). This is due to the fact that deg (e^x) — ei(x)) < 
A < ri < deg(dij(x)). Therefore, we have fci(x) = 
in(x) mod r,i(x) from (fl2l> . (fl3l >. and (II 4b . Since there 
are only t unrestricted residue errors and any other error is 
bounded in the residues, there are at most t residue errors 
with degrees greater than A. In other words, there are at least 
L — 1 — t residues dj(x) with % f 1 that are error-free or 
with bounded errors. Therefore, there are at most t errors in 
fei(x).^ 3 i(^), • • ■ ,£li{x)) to calculate k\ (x) in Step 4 of 
Algorithm I. Due to t < |_(d — 1)/2J and d < w^\ fci(x) 
can be accurately determined in Algorithm I by applying the 
residue error correction algorithm based on Propositions [2] [3] 
for the polynomial remainder code with moduli Tii(x) for 
2 < i < L, i.e., Aq(x) = k\{x). 

However, if di(x) is known as a residue with an error of 
degree greater than A, it is not guaranteed that there are at 
most — 1)/2J errors in (| 2 i(x’), • • • ,£li(x)) in Step 4 

of Algorithm I. Therefore, following the decoding algorithm in 
Section El ki (x) may not be reconstructed, and even though 
fci(x) is reconstructed, fci(x) = fci(x) may not hold. ■ 

From the above results, we have the following theorem. 

Theorem 2: Let to,(x), 1 < i < L, be L non-pairwise 
coprime moduli, M (x) be the 1cm of the moduli, and the er¬ 
roneous residue vector of a(x) with deg (a(x)) < deg (M(x)) 
be denoted as (di(x), a 2 (x), • • • ,di(x)). Denote by d the 
code distance of the polynomial remainder code with moduli 
rrii(x) for 1 < i < L. Assume that there are t<L(d-l)/2j 
unrestricted errors and an arbitrary number of bounded errors 
in the residues. Then, if the remainder error bound A satisfies 

A < T(£_ 2 |_(d- 1 )/ 2 J), (22) 

where T( t ) denotes the i-th smallest element in {ti, • • ■ ,tl}, 
each Tj for 1 < j < L is defined as Tj = inin{r,;j, for 1 < 

i 

i < L,i ^ j}, and = deg (dij(x)), we can robustly 
reconstruct a(x) as a(x), i.e., deg (a(x) — d(x)) < A, by 
following Algorithm I. 

Proof: Without loss of generality, we assume n > r 2 > 

• • • > tl- First, by taking every residue in the first 2[(d — 
1)/2J + 1 residues as a reference and following Algorithm 
L we want to calculate the corresponding folding polynomial 
fcj(x), respectively. When d i(x) for 1 < i < 2|_(ci — 1)/2J + 1 
is known as an error-free residue or a residue with a bounded 
error and the bound is A, since A < fL- 2 [(d-i)/ 2 J) < A from 
(l22b . it follows from Lemma [3 that (x) can be accurately 
determined by Algorithm /, i.e., kfx) = fci(x). Since there 
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are at most |_(d— 1 )/ 2 J residues with errors of degrees greater 
than A in the first 2 |_(d— 1 )/ 2 J +1 residues, there are at least 
[ (d — 1 )/ 2 J + 1 error-free residues or residues with bounded 
errors. Therefore, at least |_( d — 1)/2J +1 folding polynomials 
out of ki(x) for 1 < * < 2 1_(d — 1 )/ 2 J + 1 are accurately 
determined. However, when di(x) is a residue with an error 
of degree greater than A and taken as a reference, the corre¬ 
sponding folding polynomial ki(x) may not be reconstructed 
in Algorithm /, and even though ki(x) is reconstructed, it may 
not be equal to ki(x). 

Then, for each obtained k- L (x) with 1 < i < 2|_(d—1)/2J+1, 
we reconstruct a(x) as a}^(x) = ki(x)rrii(x) + di{x) for 
1 < * < 2 [(d — 1)/2J + 1. If ai(x) and tij(x) with i ^ j are 
both error-free residues or residues with bounded errors, we 
have deg (aW(s) — a(x)) < A and deg (fi^a;) — a( x)) < A, 
since ki(x) and kj(x) are accurately determined. Thus, we 
have deg (a^ (x) — (x)) < A. If di(x) is a residue with 

an error of degree greater than A and k{(x) is reconstructed, 
one can see that no matter whether fcj( x) is accurate or not, 
we will have deg (a.W (x) — a(x)) > A. This is due to the fact 
that afi(x) — a(x) = (ki(x) — ki(x))rrii(x ) + (di(x) — at(x )) 
and deg (mj(x)) > deg (di(x) — ai(x)) > A. Furthermore, 
we can easily obtain deg (aW(a;) — a^(a;)) > A when one 
of the corresponding references di(x) and dj (x) used for 
reconstruction in Algorithm I is a residue with an error of 
degree greater than A, and the other is an error-free residue or 
a residue with a bound error. 

Therefore, among the above at most 2[(d — 1)/2J + 1 
reconstructions afi(x), we can find at least [(d — 1 )/ 2 J + 1 
reconstructions such that deg(aM(:r) — a^'l (x)) < A among 
pairs of i,j with 1 < i ^ j < 2[(d — 1)/2J + 1. One 
can see that all of such reconstructions are in fact obtained 
when references dj(ir) are error-free residues or residues with 
bounded errors, and thus, any one of such reconstructions can 
be thought of as a robust reconstruction of a(x). At this point, 
we have completed the proof. ■ 

According to the above proof of Theorem [2] let us sum¬ 
marize the robust reconstruction algorithm for a given set 
of moduli {rrii(x)}fL 1 , with which the polynomial remainder 
code has the code distance d. Assume that t\>T 2 > ■ ■ ■ >tl 
and there are t < |_(^ ~~ 1 )/ 2 J unrestricted errors and an 
arbitrary number of bounded errors with the error bound A 
given by (l22l > in the residues. 

1) For every i with 1 < i < 2|_(d— 1)/2J +1, take ai{x) as 
a reference and follow Algorithm I. We want to calculate 
the corresponding ki(x) for 1 < i < 2 L(d-l)/ 2 j+l, 
respectively. Note that some kj{x) may not be recon¬ 
structed. 

2) Reconstruct a(x) as d\ l \x) = ki{x)m.i{x) + di{x) for 
each obtained ki{x). 

3) Among these obtained afi(x), we can find at least 

|_(d — 1 )/ 2 J + 1 reconstructions d^ lj ^(x) for 1 < j < M 
with > [ (d — 1 )/ 2 J + 1 such that every pair of them 
satisfy deg (x) — (x)) < A for 1 < g, q < /x, 

C / g. Then, any one of such d[ i 3 'l(a;) for 1 < j < n 
can be regarded as a robust reconstruction of a(x), i.e., 
deg (a[* 3 l(a;) — a(x)) < A. 


Remark 3: There is a related paper dealing with robustly 
reconstructing an integer from erroneous remainders [371, but 
our paper investigating the robust reconstruction problems for 
polynomials differs from E3 in several aspects as follows: 

a) The problem of robust reconstruction for integers from 
erroneous residues was considered in El and l38il . 
|[39l with different approaches. In El, a large integer 
is robustly reconstructed through constructing a new 
consistent residue vector from the erroneous residues. In 
(38], lf39l . however, all folding integers are first accurately 
determined, and then a robust reconstruction is provided 
as an average of all the reconstructions from all the 
determined folding integers. In this paper, an improved 
reconstruction algorithm for polynomials in Algorithm 
I is proposed by combining the approach in |39l with 
the error correction algorithm for polynomial remainder 
codes in 1(331 . While both of the approaches in (37l and 
(39l can be directly extended to robust reconstruction 
for polynomials in Section [Till the obtained maximal 
possible robustness bounds would be usually less than 
(ITTl ) obtained in our proposed algorithm. 

b) In (37l . a special class of residue number systems with 

non-pairwise coprime moduli was only considered, where 
moduli nrii = for 1 < i < L, dij = dji, 

and {d^, for 1 < * < L\ i < j < L} are pairwise 
coprime and greater than 1. According to Proposition 
Q] for integers, one can see that this residue code with 
these moduli m, for 1 < i < L has code distance 
2, which is unable to correct any residue errors. So, 
in order to enable single errors to be corrected, the 
legitimate range of the code must be restricted to a 
suitable subrange of [0,lcm(mi, • • • ,itil)) in li24l . and 
in El, robust reconstruction for the case of a combined 
occurrence of a single unrestricted error and an arbitrary 
number of small errors in the residues was considered 
also with the legitimate range being a suitable subrange 
of [0, lcm(mi, • • • , mr,)). Its approach is hard to deal 
with the case of multiple unrestricted errors combined 
with small errors in the residues due to a considerable 
decoding complexity. In this paper, however, we consider 
the robust reconstruction problem for polynomials from 
the perspective of polynomial remainder codes with non- 
pairwise coprime moduli. The range of the degree of a(x) 
is fixed for a general set of moduli {mi{x)}{L 1 , i.e., 
deg(a(x)) < deg (lcm (mi (x), ■■ ■ ,mi(i))). Under the 
assumption that the code distance of the polynomial re¬ 
mainder code with moduli rrii( x) for l<*<Lisd>3, 
the remainder error bound and/or the maximum possible 
number of unrestricted residue errors are obtained for 
the robustness to hold in the paper. Moreover, a well- 
established algorithm based on Theorem [2] is proposed to 
robustly reconstruct a polynomial when there are multiple 
unrestricted errors and an arbitrary number of bounded 
errors in the residues, where Algorithm I needs to be 
implemented 2[(d — 1)/2J + 1 times. 

c) Compared with El, we omit to consider the case of 
erasures in the residues in this paper due to its triviality. 
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As we know, when erasures occur in the residues, both the 
number and positions of erasures are known. Without loss 
of generality, assume that K erasures occur and the erased 
residues are aL_N+i(x), • • • , at( x). So, we just calculate 
the code distance of the polynomial remainder code with 
moduli nii(x ) for 1 < i < L — H according to Propo¬ 
sition [ 1 ] and consider to robustly reconstruct a(x) with 
deg(a(x)) < deg (lcm • • • , TOl_n(x))) from 

those available a*( x) for 1 < i < L — H in Theorem 
Q]and Theorem |2] 

Example 2: Let L = 5 and the moduli be mi(x) = ( x 3 + 

1) (x 2 —2)(x 3 +4), W 2 (x) = (x 3 + l)(x 3 — l)(x 3 +2), m^(x) = 
(at 3 — l)(x 3 + 2)(x 3 + 4), m^x) = (x 3 + l)(x 3 + 2)(at 2 — 

2) , W 5 (at) = (x 3 — l)(at 2 — 2)(x 3 + 4). Then, the lcm of 
all the moduli is M(x) = (x 3 + l)(x 3 — l)(at 3 + 2)(x 2 — 
2)(.x " 3 + 4). According to Proposition!]] the code distance of 
the polynomial remainder code with the moduli is d = 3. In 
addition, we can calculate n = 72 = 73 = 3 ,T 4 = 75 = 2 , 
and the error bound A < 3 in (l22l . Assume that there are 
one unrestricted error and an arbitrary number of bounded 
errors with degrees less than 3 affecting the residue vector of 
a polynomial a(x) with deg(a(at)) < deg(M(x)) = 14. 

1) For every i with 1 < i < 3, take di(x) as a ref¬ 
erence and follow Algorithm I. We want to calculate 
ki(x),k- 2 {x),k-z{x), respectively. Note that some fcj(x) 
may not be reconstructed. 

2) Reconstruct a(x) as a^(x) = ki(x)mi(x) + di(x) for 
each obtained fcj(x). 

3) Among aW (x) for 1 < i < 3, we can find at 
least two reconstructions at ll l(x) and d^{x) such that 
deg (x) — at* 2 l (x)) < A < 3. 

Then, a(x) is robustly reconstructed as d^(x) or a\ l2 ^(x). 

Remark 4: The robust reconstruction in Section [Til] or Sec¬ 
tion El although, may not correct all the residue errors, if 
there is an additional error correction code on the top of it, 
it may be possible to correct all these residue errors, since 
only bounded errors are left in the reconstructions due to the 
robustness. As an example for the robust reconstruction in 
Theorem |2] let the above a(x) be encoded by a product code 
in polynomial residue number system, introduced in l40l . i.e., 
given a polynomial G(x), called the generator of the product 
code, a polynomial a(x ) with deg(a(x)) < deg(M(x)) is 
legitimate in the product code of generator G(x) if a(x) = 
0 mod G(x), else it is illegitimate. Using this product code 
on the top of robust reconstruction in Theorem [2 a(x) can 
be accurately determined as a[x) = a(x) — [d(x)] G ( x ) if 
deg(G(x)) > A, where A is the remainder error bound in 
(l22l and a(x) is the robust reconstruction from Theorem 
12 This is due to the fact that deg(a(x) — d(x)) < A 
and thus a(x) and a(x) have the same folding polynomial 
(a(x) — [d(x)]< 3 (:r))/G(x) with respect to the modulus G(x). 
Furthermore, since [(t(x)] G (x) = 0, i.e., a(x) = 0 mod G(x), 
we have a(x) = a G)~[*G)\G(*) . = _ [a(x)] G(x) . 

While the above robust reconstruction has limitations in 
practice due to the type of bounded residue errors (i.e., only the 
last few coefficients of the polynomial residue are corrupted 
by errors), the theoretical result is new and may be interesting. 


In the next section, another motivation for us to study such 
bounded residue errors is shown for the improvements in 
uncorrected error probability and burst error correction in a 
data transmission. 

V. Correction of multiple unrestricted errors 

AND MULTIPLE BOUNDED ERRORS IN POLYNOMIAL 
REMAINDER CODES 

Let toi(x), ni 2 (x), ■ ■ ■ ,TOl(x) be L non-pairwise coprime 
moduli, M(x) be the lcm of the moduli, and d be the code 
distance of the polynomial remainder code with moduli to, (x) 
for 1 < i < L. As one can see in the preceding section, 
a polynomial a(x) satisfying deg(a(x)) < deg (M(x)) can 
be robustly reconstructed when there are t < L (d - 1 )/ 2 J 
unrestricted errors and an arbitrary number of bounded errors 
with the remainder error bound A given by (l22l) in the residues 
di(x) for 1 < 1 < L. Note that the reconstruction may not be 
accurate but robust and all the residues are allowed to have 
errors. Moreover, as stated in |[33ll . the polynomial remainder 
code with moduli rrii(x) for 1 < i < L can correct up to 
L(d— 1 )/ 2 J errors in the residues, i.e., o(x) can be accurately 
reconstructed when there are only t < |_(d — 1 )/ 2 J unrestricted 
errors in the residues, and any other residue is error-free. In 
this section, by making full use of the redundancy in moduli, 
we obtain a stronger residue error correction capability, that 
is, for the given set of moduli rrii(x), in addition to *<L(d- 
1 )/ 2 J unrestricted errors in the residues, some bounded residue 
errors can be corrected in the polynomial remainder code with 
moduli TOj(x) and code distance d. 

The remainder error bound r]{9) here, which depends on a 
variable 6 for 1 < 9 < L — 2[(d — 1)/2J, is given by 

V(0) < Ho), (23) 

where T(g) denotes the 0-th smallest element in {tt, • • • ,tl}, 

and Tj = minjdeg (dij(x)), for 1 < j < L, j ^ i} for 1 < 
i 

i < L. It is easy to see that the upper bound T(g) of the 
remainder error bound r](9) increases as 9 increases, i.e., the 
degrees of multiple bounded errors can be large as 9 becomes 
large. Later, we will illustrate that the larger 9 is, the smaller 
the number of correctable bounded errors is. 

Theorem 3: Let rrii(x), 1 < i < L, be L non- 

pairwise coprime polynomial moduli, d denote the code 
distance of the polynomial remainder code with moduli 
TOi(x), TO 2 (x), • • • ,TOl(x). Then, the polynomial remainder 
code can correct up to |_(d — 1 )/ 2 J unrestricted errors and 
l(L—9)/2\ — \_(d—l)/2\ bounded errors as deg(ej(x)) < tj(0) 
with the remainder error bound r](9) given by (l23l) in the 
residues. 

Proof: Without loss of generality, assume that t-\ > t 2 > 
■ ■ ■ > tl■ Similar to the proof of Theorem |2 we take every 
residue in the first L — 9+1 residues as a reference and want to 
calculate the corresponding folding polynomials by following 
Algorithm I. After that, we reconstruct a(x) as fiM(x) with 
each obtained folding polynomial as d^(x) = fc,(x)mj(x) + 
dj(x) for 1 < i < L — 9 + 1. Since rj(9) < T(g) < for 
1 <*<L — 0 + 1, ifa reference di(x) is an error-free residue 
or a residue with a bounded error and the bound is rj(9), ki(x) 
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obtained from Algorithm I is accurate according to Lemma [3] 
Since there are at most \ (d — 1)/2J unrestricted errors and 
[(L — 0)/2\ — [(d—1)/2J bounded errors in the residues, there 
are at most [(L — 9)/ 2J erroneous reconstructions in these 
(x) for 1 < i < L—6+1, and the remaining reconstructions 
are correct and equal to each other. Therefore, we can find at 
least \{L — 9)/ 2] + 1 reconstructions (x) for 1 < 3 < v 
with v > \{L — 9)1 2] + 1 such that 

a [ill (x) = a [i2] (x) = ■■■ = a w (x). (24) 

Let a(x) be equal to the majority of all the L — 9 + 1 recon¬ 
structions fiM (x) for l<i<L — 6 + 1, i.e., o(x) = a^'l (x) 
for 1 < j < v. Then, a(x) is equal to the true a(x), i.e., 
a(x) = a(x). ■ 

Remark 5: From the above proof of Theorem [3] we now 
propose our new decoding algorithm for polynomial remainder 
codes with non-pairwise coprime moduli in the following. 
Without loss of generality, assume that t\ > > • ■ • > t^. 

1) For every i with l<i<L — 6 + 1, take di(x) as a 
reference and follow Algorithm I. We want to calculate 
fcj(x) for 1 < i < L — 6 + 1, respectively. Note that some 
fcj(x) may not be reconstructed. 

2) Reconstruct a(x) as a^(x) = ki{x)rrii{x) + fii(x] for 
each obtained fc, (x). If the number of the obtained k t (x) 
is less than \(L — 9) /2] +1, the decoding algorithm fails. 

3) Among these reconstructions fiM (x) for 1 < i < L—6+1, 
if we can find at least \{L — 9)/ 2] + 1 reconstructions 
a 1^1 (x) for 1 < j < v with v > \~{L — 9)/ 2] + 1 such 
that 

a^(x) = a^(x) = ■■■ = a^(x), (25) 

let a(x) = aJ* 3 'l(x) for 1 < j < v. Otherwise, the 
decoding algorithm fails. 

With the above decoding algorithm, if there are |_(d — 1)/2J 
or fewer unrestricted errors and [(L — 9)/ 2J — |_(d — 1)/2J 
or fewer bounded errors with the remainder error bound r}{6) 
given by (1231) in the residues, a(x) can be accurately recon¬ 
structed from Theorem [2 i.e., a(x) = a(x). It is obviously 
seen that the price paid for the increased error correction 
capability is an increase in computational complexity. In the 
above decoding algorithm. Algorithm I needs to be imple¬ 
mented L — 9 + 1 times, i.e., the decoding algorithm in l33l 
(or in Section m in this paper) used to reconstruct a folding 
polynomial in Step 4 of Algorithm I needs to be implemented 
L — 9 + 1 times. 

Example 3: Let mi(x) = (x + l)(x + 2)(x + 3)(x + 4), 
TO 2 (x) = x(x + l)(x + 3)(x + 4), TO 3 (x) = x(x+l)(x + 2)(x + 

4), TO 4 (x) = x(x + 2)(x + 3)(x + 4), TO 5 (x) = x(x + l)(x + 
2)(x + 3) be L = 5 moduli in GF(5)[x]. We can easily obtain 
that the polynomial remainder code with the moduli has the 
code distance d = 4, and r, = 3 for all 1 < i < 5. Therefore, 
from Theorem 0 the polynomial remainder code can correct 
up to one unrestricted residue error and one bounded residue 
error with the remainder error bound 77 ( 1 ) < 3. If one applies 
the result in ll33ll . only one unrestricted residue error can be 
corrected. 

To see the improvements that are achieved in Theorem^ we 
consider the application in a data transmission. In the residue 


number system, a number might be communicated from the 
sender to the receiver through the transmission of its residues. 
Instead of numbers, a method for transmitting information 
based on polynomials over a Galois field is used. To simplify 
the analysis, let a sequence be a = (o[l], a[ 2 ], • • • , a[/c]), 
where a[i\ G GF(p) for all 1 < i < k and p is a prime. Denote 
by a[x) the corresponding polynomial a ( x ) = Ei=i a [*K 1 - 
Let moduli toi(x), 7712(2:), •• • , TOl(x) be L polynomials in 
GF(p)[x] such that the degree of the 1 cm M(x) of all the 
moduli is greater than k — 1. If the degree of TOj(x) is 
denoted by to* for each 1 < i < L, the corresponding 
residue aj(x) of a(x) modulo rrii{x) can be represented by 
cii(x) = dijxi- 1 for a.jj G GF (p). In place of the orig¬ 

inal block a, the residue sequences ai = (a, 1 , a* 2 , • • ■ , ai mi ) 
for 1 < i < L are transmitted in the following order: 

(ttll, , fllmi 5 ? * * * 5 tt2m2 5 * * * ; &L 1 j 5 ttLrriL ) ■ (26) 

If there is no error in the transmission, a(x) can be accurately 
recovered using the CRT for polynomials in Q. provided that 
k and the moduli TOi(x) are known. Then, a is simply formed 
from the coefficients of o(x). In practice, data can be corrupted 
during transmission. For a reliable communication, errors must 
be corrected. We herein consider two kinds of errors in the 
channel: random errors and burst errors. 

Let the channel bit error probability be 7 , the error proba¬ 
bility of a residue a* be p mi ( 7 ), and the bounded residue error 
probability of di be q mi ( 7 ; 9), where rrii is the length of the 
sequence presentation of moduli TOj(x) over GF(p). Then, 

PmA 7) = l-(l-7) mi , (27) 

q mi (T, 0) = (1- 7 ) m --^) - (1 - 7 ) m ". (28) 

In what follows, let us consider the polynomial remainder code 
in Example [2 where d = 4 and L = 5. We obtain two upper 
bounds for the uncorrected error probabilities in the decoding 
algorithm in Section [IT] and our proposed decoding algorithm, 
respectively. 

• According to Proposition 0 

1) The probability when all received residues are correct 
is 

L 

V (c) = J|(l -PmAl))- (29) 

1=1 

2) The probability when there is only one residue in error 
is 


L L 

P , ( c )='52Pm, i (.7)-'[[(l~Pm j {l))- ( 30 ) 

i=1 2=1 

oAi 

Then, from the decoding algorithm based on Proposition 
[ 2 ] in Section UU we immediately obtain an upper bound 
for its uncorrected error probability as 

Puncorrected — 1 P(^) P (t:)- (31) 

« According to Theorem [3 
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1) The probability when there are at most [(L—9) /2J = /3 
bounded errors in the residues is 


P(C ) = P(c) + ^ Qrm (7; 0) ■ n^ 1 -Prnjh)) 


i= 1 


3 =1 


L L 


+ 51 51 9m tl ^ ( 7 ; 0 ) ■ n (! — Pm, ( 7 )) 


< 1=1 <2 ><1 

L L 


3 = 1 

i#u ,*2 


EE'” E 9 m*! (7; 0) ■ ■ ■ qm ip (7; 0 ) 

<i = l < 2><1 ip>ip-i 




j=i 

,*/3 


(32) 


2) The probability when there are one error with degree 
greater than r](6) and at most [_(L — 6) /2J — 1 bounded 
errors with the remainder error bound Tj( 0 ) in the 
residues is 


L 

FR=E( 1 -( 1 - 7 r- w ) 

2—1 


L L 

+ ' 52 <lm il {r, 0 ) ■ 

*1=1 i=1 

* 1 ^* j&dl 


n^-Pmjil)) 

.3=1 


L L L 

+E E 9m 4l (7; %m* 2 (7; 0 ) ■ (1-^(7)) 

* 1=1 * 2 >*l J =1 

ilT^i *2#i j^i,il,i2 


L L L 

+ ' " + E E " ' E 9 m*! ( 7 ; (')••• 9 m i@ _! ( 7 ; 

*1=1 *2 >*l ip-l>ip-2 

* 1 #* * 2 #* 


I| ( 1 ^Pmi(7)) • 

1=1 

1 / 

(33) 


Then, from our decoding algorithm based on Theorem 0 
we immediately obtain an upper bound for its uncorrected 
error probability as 



Channel bit error probability-/ 


Fig. 1. Uncorrected error probabilities based on Proposition 
[2] and Theorem [2 simulations and theoretical upper bounds. 


We next investigate the burst error correction capability 
in polynomial remainder codes with non-pairwise coprime 
moduli. As a residue di(x) occupies to,; bits, an error in 
this residue would affect up to rrii bits. In order to express 
our question more precisely, we assume that all the moduli 
mi(x),m, 2 (x), ■ ■ ■ ,itil(x) have the same degree to. Then, 
any error burst of width not more than to + 1 in (l26l > can 
affect two residues at most. Similar to the result for polynomial 
remainder codes with pairwise coprime moduli in (29), [30 j, it 
is directly obtained that the polynomial remainder code with 
non-pairwise coprime moduli and code distance d that can 
correct up to |_(d — 1)/2J errors in the residues can correct 
up to [|_(d — 1)/2J/2J bursts of width not more than to + 1 , 
or correct one burst of width ([_(d — 1)/2J — l)m + 1. By 
Theorem [2 however, we can further improve the capability 
of burst error correction in the polynomial remainder codes 
with non-pairwise coprime moduli. Let A = [(d — 1)/2J and 
B = \_(L — 0)/2j — [(d — 1)/2J, and we have the following 
result. 




rj( 6 ) bits 


m bits 

m bits 

in bits 

m bits 

in bits 



)» + rf 0 ) bits 



r 




m bits 

m bits 

m bits 

m bits 

m bits 


Puncorrected — 1 p(t-) P (c)- (34) 


It is obvious to see that p(c ) < p(c) and p'(c) < p'(c). 
Therefore, we have 1 — p(c) — p'(c) < 1 — p{c) — p'{c). The 
performance of uncorrected error probabilities in Example [2 
for the two decoding algorithms based on Proposition [2 and 
Theorem [2 is shown in Fig. 1, where both simulations and 
the obtained upper bounds for uncorrected random errors are 
shown. 


Fig. 2. Burst error representations. 

Corollary 1: Let rrii(x), 1 < i < L, be L non-pairwise 
coprime polynomial moduli with the same degree to. Assume 
that the code distance of the polynomial remainder code with 
the moduli is d. Define by M(x) with deg (M(x)) > k — 1 
the least common multiple of all the moduli, and let rj(9 ) be 
defined in (122 . If a sequence a = (a[l],a[2], • • ■ , a[fc]) over 
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CiY(p) is encoded and sent by transmitting the coefficients 
of the residues of its corresponding polynomial a(x) modulo 
mi(x) as in (l26l >, then, based on residue error correction 
capability of polynomial remainder codes in Theorem [3] the 
following results are easily obtained: 

1) Correction of bursts of width not more than if 9) 

1.1) when A < B, it can correct up to A such bursts; 

1.2) when A > B, it can correct up to B + [(A — B)/ 2J 
such bursts. 

2) Correction of bursts of width not more than m + rj(9) 

2.1) when |_A/2J < B, it can correct up to [A/2J such 
bursts; 

2.2) when (_^4/2J > B, it can correct up to B + [(^4 — 
2B)/3J such bursts. 

Proof: From Theorem [3] the polynomial remainder code 
with moduli mfx) for 1 < i < L can correct up to A 
unrestricted errors and B bounded errors with the remainder 
error bound r](6) in the residues. Fig. 2 shows that an error 
burst of width not more than r](9) in (l26l > can, at most, give 
rise to a residue with an error and a residue with a bounded 
error. So, 1.1) and 1.2) are easily obtained. Similarly, from 
Fig. 2, an error burst of width not more than m + r](6) in 
(l26l > can, at most, give rise to two residues with errors and 
one residue with a bounded error. So, when [_^4/2j < B, it 
can only correct up to |_^4/2J such bursts. When |_^4/2J > B, 
in addition to correcting B bursts of width not more than 
m+r)(9), the remaining error correction capability can correct 
up to [(^4 — 2f?)/3j such bursts more. ■ 

Remark 6: In the previous result in 031 , a polynomial 
remainder code with code distance d can correct up to 
A = [fd 1)/2J residue errors. Based on this error correction 
capability, it can only correct up to L^4/2J bursts of width not 
more than rj(9), or correct up to |^4/3J bursts of width not 
more than m +r](9) if if 9) > 1, which are not as good as the 
above results. 


VI. Conclusion 

In this paper, we studied polynomial remainder codes with 
non-pairwise coprime moduli. We first considered the robust 
reconstruction problem from erroneous residues, namely ro¬ 
bust CRT for polynomial problem, where all residues are 
allowed to have errors, but all the errors have to be bounded. 
A sufficient condition for the robustness bound was obtained, 
and a reconstruction algorithm was also proposed in the 
paper. Then, by releasing the constraint that all residue er¬ 
rors are bounded, another robust reconstruction was proposed 
when multiple unrestricted errors and an arbitrary number 
of bounded errors have occurred in the residues. Finally, 
compared with the previous residue error correction result in 
polynomial remainder codes, interestingly, our proposed result 
shows that in addition to correcting the number of residue 
errors as in 133], some bounded residue errors can be corrected 
as well. With our proposed result in residue error correction, 
better performances in uncorrected error probability and burst 
error correction in a data transmission can be achieved. 
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